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In ABO 3 perovskites, oxygen octahedron rotations are common structural distortions that can 
promote large ferroelectricity in BiFeOs with an R3c structure [I[, but suppress ferroelectricity in 
CaTiOs with a Pbnm symmetry 0]. For many CaTiOs-like perovskites, the BiFeOs structure is a 
metastable phase. Here, we report the stabilization of the highly-polar BiFeOs-like phase of CaTiOs 
in a BaTiOs/CaTiOs superlattice grown on a SrTiOs substrate. The stabilization is realized by a 
reconstruction of oxygen octahedron rotations at the interface from the pattern of nonpolar bulk 
CaTiOs to a different pattern that is characteristic of a BiFeOs phase. The reconstruction is inter¬ 
preted through a combination of amplitude-contrast sub O.lnm high-resolution transmission electron 
microscopy and first-principles theories of the structure, energetics, and polarization of the super¬ 
lattice and its constituents. We further predict a number of new artificial ferroelectric materials 
demonstrating that nonpolar perovskites can be turned into ferroelectrics via this interface mecha¬ 
nism. Therefore, a large number of perovskites with the CaTiOs structure type, which include many 
magnetic representatives, are now good candidates as novel highly-polar multiferroic materials [^. 


INTRODUCTION 

New mechanisms to generate ferroelectricity (FE) have 
recently been the subject of active research, due to both 
fundamental interest and the technological importance 
of ferroelectrics and related materials [J]. Novel ferro¬ 
electrics have potentially higher performance for practi¬ 
cal applications, as well as potential compatibility with 
other functional properties such as magnetism, yielding 
multiferroics and other multifunctional materials [3|, l5|, |^ . 

Artificially structured perovskite superlattices offer 
rich opportunities for novel ferroelectricity IMil- Non¬ 
bulk phases for the constituent layers can be stabilized by 
the mechanical and electrical boundary conditions char¬ 
acteristic of a superlattice [mil, potentially turning 
constituents that are nonpolar in bulk form into fer¬ 
roelectrics Bill- Competing low-energy metastable 
phases can be readily found in perovskites with low tol¬ 
erance factors, promoting oxygen octahedron rotation 
(OOR) instabilities along the Brillouin-zone-boundary R- 
M line. The ground state structure in such cases is gen¬ 
erally the nonpolar orthorhombic Phnm structure. As 
a typical example, the oxygen octahedron in a CaTiOs 
(CTO) can be described by its rotation around [110] axis 
and an in-phase rotation around [001] axis (a“a“c+ in 
Glazer notation). Such a pattern of OOR favors antipolar 
behavior instead of FE Q. On the other hand, OOR with 
a different pattern can also promote large EE. As one 
famous example, in BiEeOs (BEO) with R3^c structure, 
the oxygen octahedron can be characterized by a rotation 
around [110] and an out-of-phase rotation around [001], 


yielding a fairly large polarization along [111] (a“a“a“ in 
Glazer notation). Gompared to the widespread of GTO- 
like materials, BEO-like perovskites are relatively rarely 
seen. As a result, the OOR is generally thought to sup¬ 
press EE in perovskites. 

However, for many perovskites, the BEO-like struc¬ 
ture serves as a low-energy metastable phase Q. There¬ 
fore, it would be beneficial if an artificial perovskite su¬ 
perlattice could stabilize this metastable phase for the 
entire constituent layers or in a region near interface. 
To this end, a reliable design mechanism can be derived 
only from precisely determined atomic positions in ex¬ 
periments followed by theoretical interpretations based 
on first-principles calculations. 


EXPERIMENTAL AND FIRST-PRINCIPLES 
RESULTS 

Aberration-corrected high-resolution transmission 
electron microscopy (HRTEM) is a powerful method 
for accurate visualization of oxygen octahedron distor¬ 
tions BB- Recently, it was shown that amplitude 
contrast imaging in HRTEM could be used to dis¬ 
criminate heavy and light element columns based on 
channeling contrast [l8|, allowing one to locate the 
exact interface and to visualize OOR angles in differ¬ 
ent atomic layers (see Supplementary Materials SI). 
Fig. [T](a) shows an experimental HRTEM image of 
a 4 BaTi 03 (BT 0 )/ 4 CT 0 superlattice film along the 
[110] direction of the SrTiOs (STO) substrate. This 
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FIG. 1: Experimental and simulated HRTEM images 
showing oxygen octahedral tilts in the 4BTO/4CTO 
superlatice film grown on a STO substrate along [110]. 

(a). Experimental HRTEM image using an amplitude con¬ 
trast imaging method {Cs = S/xm, Cc = l/xm, A/ = —Inm). 
BaO columns (dark dots) and CaO columns (bright dots) 
show different channeling contrast. Oxygen atomic columns 
displace differently, either upward or downward, with respect 
to the central Ti atoms in box 1, 2 and 3. 0.14nm indicates 
the spacing between Ti and O columns. Accumulated electron 
dose is 3 X 10^ electrons/nm^. (b). Simulated HRTEM image 
using the atomic positions obtained from first-principles cal¬ 
culations. The simulated image matches well with the exper¬ 
imental image except a sharper column contrast in the simu¬ 
lated image compared to the experiment one because electron 
beam induced object excitations are omitted as pointed out 
by Kisielowski et al. (see Supplementary Materials S4) M- 


image was obtained by correcting both spherical and 
chromatic aberrations to achieve amplitude contrast 
imaging conditions (Cs = 3 /xm, Cc = 1/xm). In this 
image, channeling contrast between Ca and Ba columns 
is clearly observed: atomic columns of CaO and BaO 
appear as bright and dark dots, respectively; oxygen 
and Ti columns appear as bright dots. Due to the 
interdiffusion of Ba and Ca at the interface, the intensity 
at A site varies depending on the ratio of Ca and Ba as 
discussed in detail in the supplementary material (S2). It 
is seen that BTO and CTO grow coherently on the STO 
substrate, showing the same in-plane lattice constant as 
that of STO, and elongated c-axis in the BTO layer and 
shortened c-axis in the CTO layer (see Supplementary 
Materials S3). Within the CTO layer of the superlattice 
(box 1), a strongly corrugated Ti 02 plane is observed in 
which the oxygen atoms displace upward and downward 
with respect to the central Ti atoms, corresponding to 
an OOR around [110] by 9°, comparable to that of bulk 
CTO. For Ti02 planes between two BaO planes (box 
3), alternating displacement of the oxygen atoms, and 
thus the amplitude of the OOR, is negligible, consistent 
with the fact that bulk BTO strongly resists OORs. For 
Ti02 planes between one BaO and one CaO plane (box 
2), the OOR around [110] is 3°, smaller than that in the 


interior of the CTO layers. 

For comparison, in Fig. mb) we present the simu¬ 
lated HRTEM image using the atomic positions of the 
4BTO/4CTO superlattice obtained from first-principles 
calculations. The simulated HRTEM image for the com¬ 
puted structure shows the same pattern of OOR as in 
the experiment (compare boxes 1,2, and 3 in Fig. (TJa) 
and (b)), with amplitudes of 12.5° in the CTO layer and 
5.5° at the interface. The quantitative difference in OOR 
around [110] angles from the experimental observation 
can be partly attributed to the fact that the experiments 
were performed at T= 300 K, while the ground state 
structural relaxation by density functional theory was at 
T = 0 K. In addition, in this image it is possible to discern 
the small uniform displacement of the oxygens relative to 
the Ti atoms in the Ti02 plane, which is associated with 
the spontaneous polarization of the superlattice. While 
this displacement is present in all the Ti 02 layers, it can 
be more easily identified in those belonging to the in¬ 
terior BTO layers, which do not have the corrugation 
associated with OOR. 

We use the atomic-scale information from the first- 
principles results for a detailed layer-by-layer investi¬ 
gation of the properties of the superlattice. We fo¬ 
cus our discussion on the 6BTO/6CTO superlattice, 
which allows a clearer distinction between the interface 
and interior layers; the corresponding results for the 
4BTO/4CTO superlattice are similar (see Supplemen¬ 
tary Materials S6). The computed spontaneous polar¬ 
ization is 29 iiCcm~^ along [001] and 11 jiiCcm~^ along 
[110]. The resulting layer-by-layer decomposed structural 
distortions and polarizations are shown in Fig. [2l 


TABLE I: Computed bulk properties of Rz(o), R°(o), Rxy(o), 
ADIS (A), Pooi( /xC/cm^), Piio( /xC/cm^), total polariza- 
tiou Pt( /xC/cm^), c/a ratio, aud the total euergies (eV) for 
straiued BTO aud CTO ou STO substrate modeled iu 20- 
atom supercells. Both fixed electric (E) held aud displace- 
meut (D) held bouudary couditious are cousidered, which are 
used to described the electric bouudary couditious of a per- 
ovskite iu its uature bulk or withiu au iusulatiug super lattice. 


Boundary 

Condition 

Fixed E fleld 

E = 0 V/m 

Eixed D held 

D = 29 fj.CIcm'^ 

Perovskite 

BaTiOs 

CaTiOs 

CaTiOs 

CaTiOs 

CaTiOs 

Symmetry 

P4mm 

Pbnm 

e-R3c 

Pbnm 

e-R3c 

Ri 

0 

7.7 

0 

7.6 

0 

K 

0 

0 

7.3 

0 

7.2 

Rcc^ 

0 

12.7 

11.5 

12.8 

11.7 

ADIS 

0 

0.29 

0 

0.32 

0 

Pool 

41.93 

0 

20.01 

27.81 

28.96 

Pi 10 

0 

26.40 

52.19 

24.53 

51.25 

Pt 

41.93 

26.40 

55.89 

37.08 

58.87 

c/a 

1.067 

0.964 

0.968 

0.965 

0.967 

Energy 

-162.645 

-164.158 

-164.066 

-164.106 

-164.040 


According to the dielectric slab model [1^, the struc¬ 
ture of the constituent layers of the superlattice should 
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FIG. 2: Computed local properties associated with ferroelectric and non-polar modes in 6BTO/6CTO super¬ 
lattice and schematic plots for the atomic displacements of the oxygen octahedron rotation and ferroelectric 
modes. First-principles calculations of 6BTO/6CTO superlattice showing layer-by-layer decompositions in-phase oxygen 
octahedron rotation around [001] R^, out-of-phase oxygen octahedron rotation around [001] R^, in-plane polarization Pno, out- 
of-plane polarization Pool, antipolar modes (AFE) represented by A-site cation displacements (ADIS), and oxygen octahedron 
rotation around [110] Rxy. The corresponding strained bulk values are also denoted by the open symbol in the plot, and in 
particular, the open strained bulk values in Pool are predicted by the dielectric slab model [ij. A highly polar BiFeOs-like 
interface phase in CTO is highlighted in blue. 


be closely related to those of strained bulk materials un¬ 
der the electrical boundary condition of a fixed displace¬ 
ment {D) field, imposed by the superlattice as summa¬ 
rized in Table [H Indeed, as shown in Fig. [2l the interior 
BTO layers have negligible OOR with a polarization of 32 
[iC jcvr? along the [001] direction. This is consistent with 
the structure and large polarization of strained BTO; the 
reduction from the strained bulk value of 42 fiC / crv? can 
be attributed to the electrostatic cost of polarizing the 
nonpolar CTO layer. Both bulk CTO and strained bulk 
CTO are characterized by the strong OORs due to struc¬ 
tural instabilities at the zone-boundary M and R points. 
Therefore, the interior CTO layers are dominated by Rxy 
and R\^ which are OOR around [110] and an in-phase 
OOR around [001] respectively as shown in Fig. [2l In 
addition, a large antipolar (AFE) mode develops in the 
CTO layers that can be clearly identified by the zig-zag 
movement of A-site displacement along [110] direction. 
It should be stressed that this antipolar distortion is a 
structural distortion at the X point favored by the tri- 
linear coupling due to the pattern of OOR in CTO-like 
materials. The above distortion in the interior CTO lay¬ 
ers can be clearly seen in Fig. [2l as well as in the TEM 


image in Fig. [T] (a) (see Supplementary Materials S5). 
This AFE mode was also recently pointed out to be the 
key to the suppressed FE in all CTO-like perovskites [H. 
Due to the applied tensile epitaxial strain, the interior 
CTO layers are polar along [110] direction with a mag¬ 
nitude of 26.4 iicjcrin? just like the strained CTO (2 q| . 

If the interface effect is negligible, the dielectric slab 
model can be used to predict the polarization, yielding a 
value of 22.4 \iCjew? along the [001] direction. The first- 
principles calculation gives Pool = 29.0 jiCIcw?. The 
discrepancy from the dielectric slab model suggests that 
the interface effect cannot be neglected. Such a large en¬ 
hancement of the polarization (^ 25%) is a strong indi¬ 
cation of a highly polar interface reconstruction. Indeed, 
examination of Fig. [2] reveals that the structure at the in¬ 
terface of the CTO layers differs significantly from that of 
the strained bulk CTO, with the OOR being suppressed 
at the interface of the superlattice. The AFE type dis¬ 
placement, which is driven by the trilinear coupling 0 
involving the OORs of CTO, is suppressed too. Further¬ 
more, a new structural pattern of OOR emerges at the 
interface: an OOR around the [110] axis and an out-of- 
phase OOR around the [001] axis for a TiOe sandwiched 
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between two interface CaO layers, with rotation angles 
comparable to those of the strained bulk in-phase OOR. 
This new structure pattern is exactly the same as one 
would observe for oxygen octahedron rotation in BiFeOs 
and similar perovskites with R3c symmetry. 


MICROSCOPIC MECHANISM 


Here, we propose that this change in structure at the 
interface can be interpreted as the local stabilization of a 
BFO-like structure different from that of the bulk CTO. 
As far as the topology of the oxygen octahedron rotation 
network is concerned, oxygen octahedra in both BFO and 
CTO rotate around [110]; however, BFO differentiates 
itself from CTO by its out-of-phase OOR around [001] 
instead of the in-phase counterpart in CTO. The out- 
of-phase and in-phase OOR around [001] originate from 
symmetry-nonequivalent structural instabilities at the R 
and M points respectively. This stabilization of a BFO- 
like structure in CTO layers near the interface is derived 
from the metastable polar e-R3c phase and is compat¬ 
ible with a much larger polarization than that in bulk 
CTO as shown in Table HI It has been shown that this 
phase cannot be stabilized relative to the Pbnm phase by 
epitaxial strain alone 0. However, in the super lattice, 
the suppression of the tilt angles by proximity to BTO, 
assisted by the electrical and mechanical boundary condi¬ 
tions that favor a phase with a component of polarization 
along [001], is sufficient to stabilize the structure [22|. 

To explore the stabilization of this phase more 
quantitatively, we constructed first-principles based 
models for the strained Pbnm phase (designated 
^“CTO" (Rz: R'xy, AFExy, FExy)) aud for the metastable 
e-R3c phase (E^bfo" (Rz ? R'xy, FExy, EEz)). Eacilitated 
by space group symmetry analysis, the models of both 
E‘‘cto" and E»bfo" are built through polynomial ex¬ 
pansions of the total energy from first-principles cal¬ 
culations with respect to the high-symmetry reference 
structure {PA/mmm phase) in terms of the amplitudes 
of the relevant modes. In the above, R^, R®, Rxy, 
AEExy, EExy, EEz represent the mode amplitude of 
in-phase OOR around [001], out-of-phase OOR around 
[001], OOR around [110], in-plane antipolar mode and 
in-plane and out-of-plane EE modes respectively. The 
resulting models are shown in the following for E“cto" 
and E“bfo" respectively as (see Supplementary Materi¬ 
als S7 for fitted coefficients): 


E“cto" — cliR\ + ^iRz + ^i2Rxy + ^2Rxy + ^^sAEE^y 
-j-^sAEE^y + a4EE^y -f- 64 EE^y -|- ciR^ Rxy 
+C2Ri^AFE2y + + c4R2yAFE2y 

+C 5 RxyFExy + ceAFE^yFEjjy 
+ diR^ Rxy AEExy, 

( 1 ) 



FIG. 3: Appearances of both BiFeOa-like and CaTiOs- 
like phases and the transition from one to the other as 
functions of oxygen octahedral rotations in strained 
bulk CaTiOs. Phase stabilities studied by the relative en¬ 
ergetics between BFO-like and CTO-like phase as plotted by 
AE as functions of magnitudes of oxygen octahedral rotation 
around [001] and [110] respectively. The FE mode amplitudes 
are represented by projected colors in the base as well. 


E“bfO" — OfiRz^ + + <^2Rxy + /^2Rxy T Q^sFE^ 

+ /33FE4 + a4FE2y + + 7iRf R^y 

+72Rf FE2 + 73Rf FE2y + 74R2yFE2+ 
75RxyFE2y + 76FE2FE2y + KiR°RxyFEzFExy. 

( 2 ) 

Assuming the angles of the oxygen octahedron rotations 
are tunable parameters under experimental conditions, 
we further define the functions (RL,Rxy) = 

min E“cto" (Rz? F'xy, AEExy, FExy), 

AFExy ,FExy 

and -F'-BFO'KR'z^R'xy) = 

min E«bfo"(F'z; Rxy, FExy, FEz). In order to 

FExy jFEz 

understand how the BEO-like phase can be stabilized 
relative to the CTO-like phase, we then evaluate 
AE = E’“beo" — *^“0X0" as functions of oxygen octahe¬ 
dron rotation magnitudes around [110] and [001]. The 
resulting AE is presented in Eig. [3l The total EE mode 
amplitudes are also presented by the color spectrum 
in the base plane in Eig. [H It can be seen that when 
the angles are fixed to the values of bulk CTO regions 
in the superlattice, as shown in Eig. [2] (R^ = 8.3° and 
Rxy = 12.6°), the CTO-like phase is strongly favored 
in energy. In the CTO-like phase, as shown in Eig. [2] 
and Table [H the antipolar distortion is favored over the 
EE distortion due to the large trilinear coupling term 
^ R^Rxy AEExy in Eq. [H Notably, when the amplitudes 
of OORs are reduced, the BEO-like phase becomes 
energetically more stable than the CTO-like phase as 
shown in Eig. [3l This indicates that the BEO-like phase 
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can be stabilized over the CTO-like phase when the 
OOR is reduced. When the above transition takes 
place, the OOR around [001] will change abruptly from 
in-phase rotation to out-of-phase rotation signifying 
a more drastic change in the topology of the oxygen 
octahedron network, as guided by the yellow plane at 
/S.E = 0 in Fig. m In addition to the pattern change 
of OOR, the BFO-like phase is generally found to have 
much larger polarization than that in the CTO-like 
phase as shown by the color spectrum in Fig. [H The 
much stronger FE polarization is expected, originating 
from the e-R3c phase; it can also be easily understood 
by the large four-linear coupling term ^ R°RxyFEzFExy, 
which promotes EE in both the in-plane and out-of-plane 
directions. 

This mechanism leads to the BEO-like phase that ex¬ 
ists at the interface of the BTO/CTO superlattice. As¬ 
suming the octahedra to be fairly rigid, the reduction 
of OOT is imposed by the adjoining BTO layer, which 
is strongly resistant to the OOR. A direct consequence 
of the stabilization of the BEO-like structure at the in¬ 
terface is that the polarization of the superlattice is 
greatly enhanced. Eor a particular choice of angles with 
R° = 5.7° and Rxy = 6.6° similar to those at the in¬ 
terface of the BTO/CTO superlattice, the computed po¬ 
larization of the BEO-like phase is over 54.0 [iCjcvr?^ 
which is comparable to the polarization in the bulk e- 
R3c CTO as shown in Table [D In the superlattice, the 
BEO-like phase is further favored by both the electric and 
mechanical boundary conditions imposed by the polar¬ 
ization of the BTO layer according to Eq. O Under the 
continuous displacement field along [001] direction, the 
electric boundary condition tends to polarize the CTO 
components with a larger EEz. Under the tensile strain, 
the mechanical boundary condition effectively enhances 
EExy. The larger EE^ and EExy tend to further lower 
the energy through ^ R°Rxy EEz EExy and stabilize the 
BEO-like phase. 


INVERSE DESIGN OF NEW FERROELECTRIC 
MATERIALS 

It has long been recognized that oxygen octahedron 
rotation can play different roles in perovskites promoting 
FE in BEO-like materials [l|, [23| but suppressing EE in 
CTO-like materials Q. However, the results presented 
here suggest that a transition between these two phases 
can be achieved through interface engineering in a super¬ 
lattice. In addition to improving the fundamental under¬ 
standing of these transitions, these results suggest a new 
pathway to induce EE in functional oxide materials. 

The enhanced polarization observed in the BTO/CTO 
superlattice studied here demonstrates this mechanism. 
To explore the potential of this approach, we predict a 
few more superlattices A’B 03 /A”B 03 as listed in Ta- 


TABLE II: Predicted super lattices with enhanced polariza¬ 
tions. cw?) and cw?) denote the 

computed polarizations for strained A’BOs and A”B03 re¬ 
spectively. PM{lJic/crin?) and Pca\.{lJic/cw?) are the expected 
polarizations from the dielectric slab model [3 and the com¬ 
puted polarizations from first-principles in nA’B03/nA”B03 
respectively. Polarization enhancement (Enh.) is calculated 
by (Pcai.- Pm)/Pm- Sub. denotes the proposed substrate for 
the epitaxial growth of super lattice. 


nA’B03/nA”B03 

psbulk 

^A'BOc. 

psbulk 

^A"BO.ct 

Pm 

Fcal. 

Enh. 

Sub. 

2BaTi03/2CaTi03 

41.9 

25.7 

25.5 

39.9 

56% 

SrTi03 

2BaTi03/2CdTi03 

50.8 

37.8 

34.2 

50.4 

47% 

NdGa03 

2KNb03/2NaNb03 

36.2 

~ 0 

19.8 

37.5 

89% 

DySc03 

2KNb03/2AgNb03 

36.2 

- 0 

22.0 

38.2 

73% 

DySc03 


blelH Within this category of tailored materials, one 
of the parent bulk A’B 03 is chosen to be a “CTO-like” 
perovskite with strong oxygen octahedron rotation^ re¬ 
sulting in an antipolar type (CaT 03 and CdTQqj l25| or 
an antiferroelectric type (AgNb 03 and NaNb 03 ) |26l.l27| 
ground state that is favored by the trilinear coupling 
term. The other parent bulk A”B 03 is chosen to 
have a large tolerance factor resisting oxygen octahe¬ 
dron rotation and a strong EE polarization. Similar to 
what we have already shown for the BTO/CTO exam¬ 
ple, the out-of-phase OORs around [001] are induced 
around the interface layers of A’B 03 (see Supplemen¬ 
tary Materials S8 for examples of 2 BaTi 03 / 2 CaTi 03 and 
2 KNb 03 / 2 AgNb 03 ). As a result, the overall polariza¬ 
tions of the superlattices are enhanced compared to the 
predictions from the dielectric slab model, which is equiv¬ 
alent to applying the charge continuity principle only and 
neglecting completely the possible interface reconstruc¬ 
tion. 

This approach to create new EE materials by inter¬ 
facial control can also be used to create new materials 
even where the building blocks could come only from 
nonpolar perovskites. In Table im we list a few pre¬ 
dicted IA’B 03 /IA”B 03 superlattices within this cate¬ 
gory. These interface materials also provide us a good 
opportunity to perform rigorous mode decompositions 
based on space group theory followed by a careful com¬ 
parison between the interface materials and the parent 
bulk compounds. The resulting mode decompositions 
and the local properties are also shown in Table Hill The 
A’B 03 is again a “CTO-like” perovskite with strong oxy¬ 
gen octahedron rotations. The above property is clearly 
represented by the large mode amplitudes of Qn^y and 
Qj^i which correspond to an OOR around [110] and in- 
phase OOR around [001] as shown in Table Hill Under 
such a pattern of OORs, the antipolar mode QafExy is 
favored, and EE is strongly suppressed resulting in zero 
polarization along all directions. On the other hand, 
A”B 03 is a strong “cubic” perovskite [28l-[3o| that does 
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TABLE III: Predicted new artificial ferroelectric materials based on non polar perovskites. (l)Mode decompositions (A) of Qro 
as out-of-phase OOR around [001] , Qri as in-phase OOR around [001], Qn^y as OOR around [110], Qfe^ as polar distortion 
along [001], QpExy as polar distortion along [110], and QAFExy as in-plane antipolar distortion. Their symmetries are labeled 
as , ^3 ^ respectively; (2) OOR angles (degrees) R°, R*, Ra;^ for out-of-phase OOR around [001], in-phase 

OOR around [001], and OOR around [110] (octahedron tilt) respectively; (3) The polarizations (/xc/cm^) along [001] direction 
Pool, along [110] direction Piio), and the total polarization Pt; (4) The polarization enhancement (Enh.); (5) Substrate for the 
epitaxial growth of the super lattices, KTa03(a=3.99A), MgO(a=4.2lA). (6) Space group symmetry (Sym.) of both strained 
bulk and super lattice. 


A’(A”)B03 

nA’B03/nA”B03 

Qro 

Mode Decompositions(A) 

Qr^ Quxy QfEz QFExy QAFExy 

OOR Angles(o) 

T) o ryi rj 

IXz ^xy 

Polarization(/ic/c?Ti^) 
Tool Elio Pt 

Enh. 

Substrate 

Sym. 

CdSnOs 

0 

1.22 

1.85 

0 

0 

0.79 

0 

12.0 

17.7 

0 

0 

0 

- 

KTaOs 

Pnma 

BaSnOa 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 

KTaOa 

P4/mmm 

ICdSnOs/lBaSnOs 

1.08 

0 

1.21 

0.20 

0.60 

0 

9.4 

0 

12.2 

8.9 

13.5 

16.2 

oo 

KTaOs 

Pc 

CdHfOa 

0 

1.27 

1.74 

0 

0 

0.69 

0 

12.3 

16.9 

0 

0 

0 

- 

KTaOs 

Pbnm 

BaHfOs 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 

KTaOs 

P4/mmm 

ICdHfOa/lBaHfOs 

1.07 

0 

1.17 

0.18 

0.64 

0 

9.3 

0 

12.2 

7.8 

10.8 

13.3 

oc 

KTaOs 

Pc 

CaZrOa 

0 

1.02 

1.82 

0 

0 

0.72 

0 

9.2 

16.2 

0 

0 

0 

- 

MgO 

Pcmn 

BaZrOs 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 

MgO 

P4/mmm 

ICaZrOa / IBaZrOa 

0.79 

0 

1.22 

0.27 

0.71 

0 

6.6 

0 

10.8 

11.1 

27.6 

29.7 

oo 

MgO 

Pc 

ICaZrOs / IBaZrOs 

0.90 

0 

1.16 

0.30 

0.67 

0 

8.1 

0 

10.6 

13.5 

22.5 

26.2 

oo 

Relaxed 

Pc 


not display structural distortions associated with either 
OORs or polarization at its ground state as shown in 
Table Hm 


Strikingly, when A’BOs and A”B 03 form a 
1 A’B 03 / 1 A”B 03 superlattice, the resulting struc¬ 
tural distortions are significantly different from their 
parent bulks. The differences come not only from the 
amplitudes of the modes but also from the symmetries 
associated with these modes. In Table mu the OORs 
around [110] Qn^y are preserved in all these superlattices 
but, with largely reduced mode amplitudes compared 
with those in bulk A’B 03 . In contrast, the in-phase 
OOR around [001] completely disappears and is 

replaced by a large mode amplitude Qro associated with 
an out-of-phase OOR around the same axis in all the 
predicted new materials. As we have seen repeatedly in 
the previous discussions, such a new pattern of OOR 
signifies the stabilization of a “BFO-like” structure in 
all these artificial materials. Accordingly, large polariza¬ 
tions develop along both [001] and [110] directions with 
the generated total polarization vector roughly along the 
[111] direction due to the broken symmetry in the e-R3c 
phase. It can be noted that the polarization of BiFe 03 
is exactly along [111] direction in the R3c symmetry. At 
the same time, the antipolar mode QAFExy completely 
eliminated. Here, we want to stress that none of the 
component perovskites in the predicted superlattices is 
polar either in its natural bulk or in its strained bulk 
formats! 


OUTLOOK 

Currently, there are two widely adopted interface ap¬ 
proaches to induce EE in oxide super lattices, namely 
tricolor 0 and hybrid improper methods [3l|. An ar¬ 
tificially induced broken inversion symmetry lies at the 
heart of both the above two methods. In the former, the 
broken inversion symmetry along the out-of-plane direc¬ 
tion is introduced by the number of species in the su¬ 
perlattice; while in the latter, the broken inversion sym¬ 
metry along the in-plane direction is facilitated by the 
differences in the antipolar modes of the two perovskite 
materials across the interface. However, it should be 
noted that the interface approach discussed here is a new 
route that is conceptually different from the above. In¬ 
stead of introducing artificial inversion symmetry break¬ 
ing, the ferroelectric polarization is stabilized by favor¬ 
ing a “BFO-like” structure which is a metastable phase 
for many perovskite materials. Due to the nature of 
the energy term that stabilizes the “BFO-like” struc¬ 
ture (^ R°RxyFEzFExy) , it is expected the switching of 
EE does not necessarily require switching the directions 
of oxygen octahedron rotations, which usually requires 
much larger energy as is implied in the hybrid improper 
mechanism. Indeed, the EE polarization switching has 
already been successfully demonstrated in 2BTO/2CTO 
by Lee’s group [32[. Based on nudged elastic band (tran¬ 
sition state) theory mill and single domain assump¬ 
tion, the energy barrier in switching EE in 2BTO/2CTO 
(154 meV) is found to be close to that of the predicted 
materials lCdSn 03 /lBaSn 03 (119 meV) both of which 
are modeled in 40-atom supercells. 

In conclusion, by combining HRTEM experimental and 
first-principles approaches, we introduced a comprehen- 
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sive interface design method to stabilize a highly po¬ 
lar “BFO-like” metastable phase in perovskite materials. 
Both the electric and mechanical boundary conditions 
are taken into account as well. This scheme introduces 
a conceptually novel way to design artificial FE materi¬ 
als. By predicting some new materials, we demonstrate 
this approach of exploring novel functional materials. 
For example, if the FE could be recovered in orthogonal 
REe 03 (R= Y,Gd, Tb Dy,Ho,Er,Tm,Yb, Lu) (aMg by 
this approach, the synthesis of a new family of room tem¬ 
perature multiferroic materials could be achieved. Eur- 
thermore, the result of our current work indicates that, 
through an interface design mechanism, short-period su¬ 
perlattices can have stronger EE than longer ones. This is 
promising for modern device applications based on ultra- 
thin films. 
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